Effects of anisotropic spin-exchange interactions in spin ladders 



R. Citro a and E. Orignac 6 
a Dipartimento di Fisica "E.R. Caianiello" , 
University of Salerno and Unita INFM of Salerno, Baronissi (Sa), Italy 
b Laboratoire de Physique Theorique, CNRS UMR 8549, 
Ecole Normale Superieure, 24 Rue Lhomond 75231 Paris Cedex 05, France 

(February 1, 2008) 

We investigate the effects of the Dzialoshinskii-Moriya (DM) and Kaplan-Shekhtman-Entin- 
Wohlman-Aharony (KSEA) interactions on various thermodynamic and magnetic properties of a 
spin 1/2 ladder. Using the Majorana fermion representation, we derive the spectrum of low energy 
excitations for a pure DM interaction and in presence of a superimposed KSEA interaction. We 
calculate the various correlation functions for both cases and discuss how they are modified with 
respect to the case of an isotropic ladder. We also discuss the electron spin resonance (ESR) 
spectrum of the system and show that it is strongly influenced by the orientation of the magnetic 
field with respect to the Dzialoshinskii-Moriya vector. Implications of our calculations for NMR and 
ESR experiments on ladder systems are discussed. 



I. INTRODUCTION 



In recent years low dimensional quantum spin systems have attracted great interest due to their fascinating prop- 
erties originating from low dimensionality and quantum fluctuations. Quantum effects are particularly relevant in 
S = 1/2 systems. The quasi-one-dimensional antiferromagnetic S = 1/2 ladders have been extensively jjwestigated 
both theoretically and experimentally^. Several experimental realizations of such systems are availableau. In par- 
ticular |-Sr n _iCu n -i-i02n and Sri4Cu24C>4i seem to be quite well described by. p. simple isotropic Heisenberg ladder 
modelso. However, countless experiments on higher dimensional compoundsETcl indicate that in many systems the 
isotropic exchange interaction is not sufficient to describe magnetic properties such as weak ferromagnetism. From 
the theoretical point of view, more than thirty years ago DzialoshinskiiEl pointed out that a term (allowed by sym- 
metry in non centrosymmetric structures) of the form D • (Si x S2), where D is the constant Dzialoshinskii vector 
and Si, 2 are the sublattice magnetizations, favors a canted spin arrangement over the antiferromagnetic one and a 
weak ferromagnetic moment. This term was later derived by Moriyaliil from ajiiicroscopic Hubbard- type Hamiltonian 
by including spin-orbit coupling in the Anderson superexchange calculationEj. Actually, the calculation of Moriya 
showed that besides the Dzialoshinskii-Moriya (DM) term, a symmetric anisotropy term of the form S ■ A ■ S was 
also obtained, but was. assumed by Moriya to be negligible compared to the antisymmetric one. However, it has 
been realized recentlytjTlj that this assumption is incorrect if the underlying microscopic model from which spin-spin 
interactions are derived has an SU(2) spin symmetry. In that case, the effective spin-spin interactions have the same 
SU(2) invariance as the microscopic model. As the DM term alone breaks SU(2) symmetry, the symmetric term has 
to be such that it compensates exactly the spin anisotropy induced by the DM term. For this class of systems, the 
symmetric anisotropy term is called tlre-Kaplan-Shekhtman-Entin-Wohlman-Aharony (KSEA) and it is tuned in such 
way that SU(2) symmetry is recoveredEjEj. A derivation of the DM and .the KSEA term starting from a microscopic 
SU(2) invariant Hubbard model with spin-orbit coupling can be found ir£j. 

The main objective of this paper is to study the influence of the anisotropic (DM and KSEA) spin-exchange 
interactions in a spin 1/2 two leg ladder. We will consider a uniform Dzialoshinskii-Moriya interaction in contrast 
with the staggered interaction discussed in [r| We mainly concentrate on a DM interaction along the rung of the 
ladder, but we discuss also briefly the case of a uniform DM interaction along the chains. Our motivation for restricting 
mainly to the case of interactions along the rungs is the situation in CaCu203 in which the Cu-O-Cu bond angle, in 
the ladder rungs is 123° indicating that in addition to the Heisenberg coupling a Dzialoshinskii-Moriya interaction^] is 
also present in the rungs of the ladder. Note that a model of coupled chains with Dzialoshinskii-Moriya interactions, 
motivated by experiments on CsCuCLi has also been considered in Ref. [ll| However, in Ref. [l8[ the authors are 
considering an two dimensional array of coupled chains using bosonization and RPA rather than an isolated ladder. 
The paper is organized as follows. We first outline the bosonization proceduE£,of a weakly coupled two leg spin ladder 
in presence of DM and KSEA interactions along the rungs. The derivatiorO of the low energy theory is now quite 
standard, and the reader already familiar with bosonization techniques should read it only to get acquainted with 
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our notations. Expressions for the staggered correlation functions can then be obtained using a mapping onto four 
non-critical Ising modelscJ. Then we study all uniform magnetic correlation functions and calculate the temperature 
dependence of the spin susceptibility and the NMR relaxation rate. Finally we discuss the very interesting effects of 
the anisotropic spin interactions on the electron spin resonance (ESR) spectrum, where the DM interaction triggers 
a phase with gapped but incommensurate correlation functions with applied magnetic field perpendicular to the DM 
vector, while the effect of a superimposed KSEA interaction is the closure of the gap with a restoration of the SU (2) 
symmetry. 



II. GENERAL CONSIDERATIONS AND BOSONIZATION 



A. General considerations 



The Hamiltonian of an antiferromagnetic two leg ladder with both Dzialoshinskii-Moriya (DM) and Kaplan- 
Shekhtman-Entin-Wohlman-Aharony (KSEA) interactions along the rungs is: 

H = J,, 2(S M Si )<+ i + S 2 ,iS 2 , m ) + Jj_ Si,iS 2 ,i + D • £(S M x S 2)i ) + s i,i*S*,i i 1 ) 

i i i i 

where D is the uniform Dzialoshinskii vector associated with the bond between the spins on different chains. The 
last term in the Hamiltonian (Q) describes the KSEA interaction, where A is the anisotropy tensor. Generally, in a 
system having full SU(2) symmetry, D and A are not independent; they are related to each other in such a way that 
a rotation of the spin axes, maps the single-bond Hamiltonian onto an isotropic one with no preferred direction of the 
staggered magnetizationli3. If we use the quantization axis z for the spins such that D = Dz, the anisotropy tensor 
reduces to A zz = A and all other components become zero. The Hamiltonian (Q) reads: 

H = J2J\\ (SmS m+1 + S 2 ,S 2 , l+1 ) + ^±^s+ s- + ^_^s M S+ + (Jj_ + A)Sl 1 Sl 2 (2) 

i 

This Hamiltonian can be simplified by the gauge transformation!!!!: 

where 2a = arctan(Z)/ Jj_), and brought to the form: 

H = J2 J||(S M S M+ i + S 2)i S 2 , i+1 ) + J±(Sf tl Sf t2 + %%) + (Jj_ + A)Sl x Sl a (4) 

i 

where Jj_ = \J J\ + £> 2 sign( J±). The KSEA interaction is given by A = J\ + D 2 sign(Jj_) — Jl leading to the 
restoration of SU(2) symmetry. For D <C Jj_, one recovers the approximate expression A ~ D 2 /(2J±) (see Eq. (2.19) 
of Ref. |l^). Let us mention that a Dzialoshinskii-Moriya interaction along the chain direction: 

^DM = D|| ' X ^» i + 1 .P (^) 

n,p 

can also be removed by another gauge transformation: 

$n,p — e%Sn Sn,pi (6) 

where 5 = arctan(D /parallel / J) . This gauge transformation has no effect on the interchain coupling. However, in the 
absence of KSEA interaction, it can induce an anisotropic interaction along the chain direction. Further, this term 
induces incommensuration in the correlation functions. Let us discuss briefly the limit = 0. With Dzialoshinskii- 
Moriya interactions only, and for classical spins, the spins Si, 2 would lie in the ground state forming an angle of 7r with 
each other. Thus, the real spins would make an angle of ir — 2a with each other. In that case, a weak ferromagnetic 
moment proportional to sin a would appear in the ground state. In the case of quantum spins 1/2, we can rewrite 
the interchain interaction up to an unimportant constant as: 

^interchain = ~2~('^ 1 + ^ 2 ) 2 ^ ^ 
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As a result, the total spin S — Si + §2 and the spin projection along the z-axis S z — S z + 5| are good quantum 
numbers. In the case of an antiferromagnetic interaction, the ground state is the singlet state with \S = 0) and the 
excited states are \S = 1;S Z = ±1) with energy E(S = 1,S Z = ±1) = J± % J± , and \S = 1;S Z = 0) with energy 
E(S = 1, S z — 0) = J±. We notice that the state with S z = is the highest excited state in this case. 

For a ferromagnetic interaction, the state \S = 1; S z — 0) becomes the state of lowest energy. The excited states are 
then \S = 1; S z = ±1) and 1/5 = 0), the latter one being the highest excited state. Clearly, in the case of ferromagnetic 
interactions, the system becomes equivalent to a XXZ spin 1 chain with easy plane interaction — K(S Z ) 2 . This could 
lead to the observation of an XY2 phasdrJ. 



B. Bosonization treatment 



For weak interchain coupling J± and D,A <C J\\, the Hamiltonian (Q) can be bosonized following Ref. [ll| 

For the moment let us consider A = 0. For J± = D = 0, the low energy properties of the Hamiltonian (Q) are 
described in terms of the boson operators <p a [a. = 1, 2) and their conjugate momenta ttH — d x 9 a . The Hamiltonian 
of chain a is: 



i.d x (t> a f 



(8) 



where u = ? Jna is the "spinon" velocity and a is the lattice spacing. Let us now turn a weak non zero J±, D in the 
Hamiltonian (||). By using the continuum limit of spin operatorsEj: 



5 n 



S^(x = na) 



b^(x — na) = — 



-A sin V2<p a 



e l " A + cos V2(j) a 



(9) 



we can rewrite ([!]) as: 

H = H s + H a 
dx 
2^ 
dx 
2^ 



H a = 



u [(7rn s ) 2 + (d x cp s ) 2 ] - 
u[{nii a f + (d x 4> a ) 2 ] + 



J ± X 2 



2ir 2 a 



J 1 A 2 



2ir 2 a 



dx cos 2(f> s 



dx cos 2<f) a 



J 1 A 2 



dx cos 29 a , 



where we have introduced the symmetric and antisymmetric combinations: 



(10) 



(11) 



V2 ' Ta V2 ' 

and similar combinations for n and A ~ 1 is the expectation value of the charge operatorEl. Note that 

the Hamiltonian H a , H s can also be obtained by bosonizing directly the Hamiltonian (|l|) and then performing a 
shift of the field 9 a which is equivalent to the gauge transformation Eq. (||). In the presence of a Dzialoshinskii- 
Moriya interaction along the chains, the effects are two-fold. First, we have to make a shift 9 P — > 9 V + S^, i.e. 
9 S + \/2/3- ■ Since in the Hamiltonian ([l(]), only the term cos20 s is present, such shift does not affect the spectrum 
of the ladder. However, it induces incommensuration. The spin-spin correlation function (S + (q)S~ (—q)) will present 
divergences for q = S/a and q = ±ir/a + S/a rather than respectively q = and q = n/a. The second effect of 
the longitudinal Dzialoshinskii-Moriya interaction is to induce an anisotropic interaction 5J z ^2 ip S z p S z +lp . Such 
interaction gives marginally relevant interactions upon bosonization. On the other hand, the interchain coupling gives 
relevant interactions of dimension 1. Thus, for not too weak interchain coupling, the extra anisotropic interaction 
will not affect the physical properties of the ladder, and the only effect of the longitudinal Dzialoshinskii-Moriya 
interaction will be a shift of the copelation functions. The Hamiltonians H a , H s can be further rewritten in terms of 
non-interacting Majorana fermionsE3: 







nil 



sin(6> s 



cos 



(9a - 4>a) 



£-R,4 — 



sin 



(9a - <j>a) 



(12) 
(13) 
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and similar definitions for the £l's with —(f> SM — * (f>s,a- Reexpressed in terms of Majorana fermions the Hamiltonian 
reads: 
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a=l J Z 

The spectrum is thus composed of a Majorana doublet (y — L,R), with mass mi,2 = m = X 2 J±/2n and 

two singlets £jj> ^ °f respective masses 7713 = A 2 (^ — ^-), 1TI4 = — A 2 (^ + ^-). Without the DM interactionEl thr; 
spectrum would be formed of a triplet Q, a = 1, 2, 3 with mass m and a singlet £^ with a larger modulus mass 3|m|£2l. 
The loss of SU{2) symmetry in the presence of a DM interaction thus translates into the lifting of the degeneracy 
of the triplet. Let us note that the Majorana doublet describes the excitations of spin S z = ±1, whereas the singlet 
corresponds to the excitations of spin S z = 0. We notice that the hierarchy of energy scales as well as the degeneracy 
of excitations with S z — ±1 that we had obtained in the strong rung coupling limit are preserved in the limit of a 
small rung coupling. In the presence of a KSEA interaction, the bosonization treatment of ([I]) is identical to the one 
in [l9| We find that mx^.z = \ 2 J±/2ir, thus recovering the triplet of the isotropic ladder. 

From the Hamiltonian Eq. (H), the specific heat of the system at low temperature is obtained in the form: 



/mi\ 3 / 2 mi /m 2 \ 3 / 2 m 2 _na , /"rn 3 \ 3 / 2 m 3 _ s 



2nu \ T / \Z2ttu \ T / \/2itu 

If the KSEA interaction is different from zero, the SU(2) symmetry is restored and the specific heat becomes, C v ~ 

( A // \ 3/2 A // a" 
■^J -^==-e~~, where A" is the triplet mass \ 2 Jj_/2tt. The difference between the regular ladder and the ladder 

with only DM interaction is the contribution from the singlet excitation with S z = 0. The specific heat as a function 
of temperature is plotted on figure [jj A remarkable feature is that at low temperature, the specific heat of the ladder 
with DM interaction is lower while it is higher at high temperature. 



III. STAGGERED CORRELATION FUNCTIONS 



To calculate the two-point correlation functions of the staggered magnetization, we use the mapping onto a pair 
of non-critical 2d Ising model introduced in Refs. This mapping permits to express the staggered correlation 

functions -terms of the correlation functions of the order and disorder parameters of the Ising model away from 
criticalityETI^J. To precise the correspondence, we will use the convention of Ref. |l9| that J±/t > 0, where t is the 
reduced temperature of the corresponding noncritical Ising model, t — (T — T c )/T c . With this convention, J± > 
corresponds to the disordered phase and we have the following bosonization formulas for order and disorder parameters 
oi and fa (i = 1, . . . ,4): 



COs((j) s ) 


= Ml|«2, 


sm(cj) s ) 


= Ol<72) 


cos(0 s ) 


= CTl^ 2 , 


sin(6' s ) 


= Ail02, 


cos(0 a ) 


= M3 °4 ; 


sm((/) s ) 


= 03^4, 


cos(6» a ) 


= C3CT4, 


sin(0 a ) 


= 



(16) 



Using formulas (|1E|) along with the gauge transformation (|3j) , we obtain the components of the total M and relative 
m staggered magnetization as: 

M x ~ (cos a <J\[iiOzo^ — sin a cfi^HzHa) Wi ~ (sin a /Xi<72<73cr4 + cos a ^i<J2^Hi) 
My ~ (cos a [i\02(Jz<Jt± — sin a /Lti<T2/i3/i4) m y ~ (sin a <j\iii(Jz<j± + cos a <7i/i2/i3/X4) 
M z ~ 0-102//3CT4 m z ~ ^1^20-3^4. (17) 

Comparing these expressions with those obtained by Shelton et a/El, we see that the DM interaction mixes the total 
and relative magnetization of the isotropic case in the x — y plane, while the z com pon ents are unchanged. This 
is a consequence of the easy-plane effect induced by the DM interaction. From Eqs.(|17|) the two-point correlation 
functions perpendicular, and along the DM vector, are expressed in term of correlation functions of the order and 
disorder parametersOt3: 
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(T T M x (x,r)M x (0,Q)) = G a ( — )G IM (— )[cos 2 a G ff (^)G CT (^)+sin 2 aG M (^)G p (^)] 



u u u 

,mr , „ ,mr, r . „ i m 3 r \^, , m 4 r ^ 



(T T m x (x, t )m x (0,0)) - G CT ( — )G M ( — )[sin 2 a G CT (^)G CT (- 



(T T M z (x,r)M z (0,0)) = [G CT ( — )] 2 G A1 (^)G CT (^) 



It 

(T r TO x (x,r)m 2 (0,0)) = [G /1 ( )\ G a { 



u 



a 



u 



cos 2 aG, 



u 

, m 3 r 



)G A 



u 



(18) 



where r = x 2 + (tir) 2 and G a {vnirju) and G^imirju) are the two-point correlation functions given by Eqs. (38) 
and (39) of Ref. |l9|. In order to determine explicitly the correlations length in the direction orthogonal and along the 
DM vector, we use the asymptotic expansions of the functional Kq(v) and K%(v) for v — > oo: 



K Q (v) = v^/(2^)e-"(l - l/(8v) + o(l/v)) 
Ki(v) = ^n/(2v)e- v (l + 3/(8v) + o(l/v) . . .). 

Substituting these expressions in Eq.(O) we get, after some cancellations: 



(T t M x (x,t)M x (0,0)) 
(T T m x (x,T)m x (0,0)) 

(T t M z (x,t)M z (0,0)) 
(T T m z (x 1 T)m z (0,0)} 



Ato sin 2 a r/u 



y / 2i:mr/u 
Ai, cos 2 a 

t 

y / 2nmr/u 



-mr/u 



At 



2ir(m 2 m4) 1 / 3 r 

— mr/u 



3/2 



(19) 

(20) 
(21) 

(22) 

(23) 
(24) 



We see that with only Dzialoshinskii-Moriya interactions, the in-plane and out of plane correlation functions have 
a different correlation length in contrast to the case of the regular ladder. This difference of correlation lengths 
should be observable in neutron scattering experiments. The ratio of the correlation le ngths should give access to 
the intensity of the Dzialoshinskii-Moriya interactions, since it is equal to m^/m — 2y/l + (D/J±) 2 — 1. Another 
difference with the regular ladder is that (M X M X ) decays as slowly as (m x m x ) with DM interactions whereas its decay 
is much faster without them. The ratio (m x m x ) / (M X M X ) provides a second measurement of D/J±. When we take 
into account the KSEA interaction, the SU(2) symmetry is restored and correlation functions have the same decay 
with r in all directions with £ = ( A 2 Jx Nevertheless, in presence of a KSEA interaction the decay of (M X M X ) 
remains not faster than the decay of (m x m x ), and the ratio of these correlation functions gives access to the strength 
of the Dzialoshinskii-Moriya interaction. Using the above expressions ( |l8|) we can easily calculate the imaginary part 
of the dynamical transverse spin susceptibility for very small momentum. We find: 



Imx XX (qi|,q_L = 0,w) 



Imx xx (q||,q_L = tt,u/) = 



A^ sin 2 a7r 2 u 
i^\J (uq) 2 + m 2 
Ai, cos 2 «7r 2 u 



Imx zz (q||,q_L = 7T,w) 



V M 2 + 

A^7T 2 U 

*V( u <l) 2 + 



6(u 


-VM 2 - 


- m 2 ) — 5 {ui - 


- y (uq) 2 + m 2 ) 


S(uj 


- V(uq) 2 - 


- m 2 ) — 8(u> - 


f- V(uq) 2 +m 2 ) 


5{uj 


- y(uq) 2 - 


hm 2 ) -5{oj- 


- \J (uq) 2 + m 2 ) 



(25) 

(26) 

(27) 
(28) 



In presence of a KSEA interaction the asymmetry disappears and the same optical magnon peak is obtained in 



both directions at uj = y {uq) 2 + (A 2 J±/2ir) 2 . Thus the analysis reveals that there is net way of determining the DM 
interaction and the interplay between DM and KSEA interactions into experiments. 
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IV. MAGNETIC SUSCEPTIBILITY AND NMR RELAXATION RATE 



A. General formalism 

To determine the temperature dependence of the -magnetic susceptibility and of the NMR relaxation rate 1 /T\ , we 
only need the slowly varying part of the spin density! 3 !]. That component of the spin density can be expressed in terms 
of Majorana fermions. The corresponding correlators can the be obtained by the approach of Ref. [30| The following 
relations hold for the sum and the difference of the spin density currents of chains 1 and 2 along the spatial directions 
(l,2,3)(see also Appendix [A|): 

(J v i + Jul) 1 = -«(cosa £,(,1 - sin a (J„i - J v2 ) 1 = -i(cosa + sin a (1&) 

(Jul + Jv2) 2 = -i{cosa tltl + sin a CO [Jyi - J^f = -i(sina + cos a ££) 

( Jul + J„2? = -itltl (J„i-Jv2?=i££, (29) 

where v — L, R. From the previous expressions we evaluate the slowly varying component of the Matsubara spin-spin 
correlation functions which are useful for neutron scattering experiments and NMR relaxation rate: 

X ah {q,iuJn)= dTdxe l( -^ T - qx) (T T J a (x,T)J b (0,Q)) T (a,b = 1,2,3). (30) 



The finite-temperature correlations are obtained by the analytic continuation iu> n — > lu + i0+ . Note that expression 
( |30| ) corresponds to a tensor susceptibility, and J a ^ is the total spin current. 



j»w= Yl (jT + jT)- (3i) 



v=R,L 

The explicit expression for the correlation functions of the uniform currents in real space is: 
(T T j\x,r) J 1 (0,0)) =- [cosM^a^)^(0,0))+sin^^^^ 



V=L,R 
H = L,R 



(T t J 2 (x,t)J 2 (0,0)) = (T t J 1 (x,t)J 1 (0,0)) 

(T t J 3 (x,t)J 3 (0,0)) = - J2 ( T t££&t)ZXM), (32) 



v = L,R 

i±=L.n 



the other contributions being zero. 

From (||^) we can evaluate the uniform susceptibility by applying Wick's .theorem. The corresponding expressions 
can be written in compact form introducing the thermal Green's functional for left- and right-moving triplet and 
singlet Majorana fermions: 



G RR (k,iw n ) = G1 L (-k,iuj n ) 



iujr, + uk 



G a RL (k,iu n ) = G? R (k,iu n ) = - + ™£ + , (33) 

where a stands for (1,2) in the case of the doublet excitation or for (3,4) in case of a singlet excitation, and the 
" polarization" function: 

T^(q, iu n ) = -/T 1 E / § t G ^.( fc ' G ^(9 ~ *, iw„ - iJ n ) + G a RL (k, iu>' n )G RL (q ~ k, zc„ - «/ n ) 

U)' 

+G a LR (k,iJ n )Gl R (q - k,iu> n - icj' n ) + Gl L (k,ioj' n )G{ L (q - k,iu n - icj' n )}. (34) 

We find: 

X 11 (q 1 ioJ n ) = [cos 2 ar 23 (g,w„) - sin 2 aT 2i (q, uj n )\ , 
X 22 (iw„) = x X1 (iu n ), 

X 33 (iLo n )=T 12 (q,cj n ) (35) 
X 12 (iu n ) =x 21 (^«) = 0. (36) 
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B. Calculation of the magnetic susceptibility 



To obtain the static susceptibility, we need to take the limit u> — > 0, and the k — > in x(A;, uS), The evaluation of 
the Matsubara frequency sum for each term of type ( |34| ) , in the limit q — > 0, leads to the following expression in the 
static limit. If m a — mp: 



1 



T 



dk 



°) = 7F ^ sech ( 



2tt 



(37) 



where t a{j3 ) k 



l 2 k 2 



m 



For m a ^ m^, a more complicated expression is obtained (see Appendix M). In this 



expression, a new contribution can be isolated 

1^(5-0,^ = 0) = 
This contribution does not vanish for T 



dk 
2^ 



1 



(uk) 2 + m a mp \ 1 — n F (e a (k)) — n F (ep(k)) 



e a {k)ei3{k) 
0, but gives: 



1 



rS /3 (g-»0,w = 0) T =o = — 



m a mp 



e a (k) + ep(k) 



ln[^ 



771 



/3 



(38) 



(39) 



We obtain the final expressions of the static and uniform components of spin susceptibility for T — > as: 



X U (0)=X 22 (0) = - 
iru 



cos 2 a 



(m§ 



2^ ln ( — ) 



. 2 

sin a 



(ml 



2^ ln 

m z ) m 



0{T 



-l/2 e -mm{m«}/T\ 



x dd (0) 



_ T -l/2 e -^_ 



(40) 



The result shows that the in-plane static susceptibility does not vanish for T — 0, whereas it decays exponentially, 
T~ 1 / 2 e~ x J -i-/( 27rT ) along the DM vector. A similar result holds ip_the case of the anisotropic spin 1 chain where the 
in-plane magnetic susceptibility is also finite at zero magnetic fieldEa. Restoration of the SU(2) symmetry by the KSEA 

interaction gives exponentially decaying susceptibilities in all directions: x n (0) = X 22 (0) = X 33 (0) — T -1 / 2 e _ ~T~ , 
where A" = X 2 J±/2%, similarly to an XXX ladder. This is shown in Fig.|| where a remarkable feature is that at low 
T the spin susceptibility is lower in presence of a DM plus KSEA interaction, with a larger spectral gap opening. 



C. NMR relaxation rate 



In this section, we calculate the NMR longitudinal relaxation rate T-y . By definition, T x 1 is related to the 
imaginary part of the dynamical susceptibility in the longitudinal direction: 



(^VT^y^M, (41) 



k 



The susceptibility to) can be decomposed into two contributions, one form the uniform part of the magnetization, 
Xo(k,uj) and one from the staggered component XttQZjW). It is well known that in the case of a ladder, the staggered 
component Xm of the spin density makes at low temperatures a negligible contribution to T^ 1 compared with the 
uniform partEj xo> so we can focus on this latter contribution which is easily obtained from the method of Ref. |3(]. 
In the following, we will denote xo(fc>k>) simply by x to simplify notations. We have: 

f 3 ■ 

X a (k,iu; n )= e^ T {J a (k,T)J a (0,0)}, (42) 



and J a (a = 1,2) is the total spin current in chain 1 or 2. It can be derived by Eq.s( A2)-( A4)-( A5) in the Appendix 
[A|. Using the polarization function, we get the following expression for the longitudinal susceptibility: 

X a {q, ioJ n ) = [- cos 2 a(r 23 (q, «j n ) + r u {q 7 iui n ) + T 3 \q, iu n ) + T 2 \q, iuo n )) 
- sin 2 a(r 24 ((7, kj n ) + T 13 (g, iu n ) + T 14 (q, iu n ) + T 23 (q, iw n )) 

-(r 12 (q,iu n )+T 34 (q,iu n ))]. (43) 
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Taking into account that the fermionic branches 1,2 have the same spectrum, we can reduce the above expression to 
a simpler form: 



X a {qMn) = -[2(r 23 (q,iuj n )+r 14 (q,iLj n )) + (T 12 (q, ioj n ) + T 3i (q, iv n ))]. 



(44) 



The [detailed calculations are reported in the Appendix IB]. Here, we only quote the final results. For m a ^ mp, we 
haveka 



TV lim 



Imr a ^(q,w) 



max(m a ,m^) 



47T 



i a mp 



sech 2 (^). 



(45) 



For m a — mp, the NMR relaxation rate is divergently for lu = 0. For ui ~ 0, the following expression is obtained: 

^m'm.^j^:^ (46) 
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where 7 is Euler's constantS. 

From (f45|) and (f4§), the final expression for the longitudinal relaxation rate reads: 



1 



de 



m 3 4tt^2 
e- m / T ln 



! -\/e 2 — to 2 . 





r°° de 


e 2 + mm4 




f im 4 4^ 


Ve 2 — m 2 e 2 — m| 



ni3ni4 



sech 2 ( — 
2 



7T 



(47) 



The first two terms give the in-plane singlet-triplet contributions, while the last two terms are the contributions along 
the DM z-axes. Compared with the results in Ref. ^ the in-plane singlet-triplet contributions are different. For a 
DM plus KSEA interaction these contributions become equal, leading to same result of Ref. ||(] with J±—*J±, where 
an explicit dependence on DM interaction is still present. 

V. CORRELATION FUNCTIONS IN APPLIED MAGNETIC FIELD AND ESR SPECTRA 

The electron-spin resonance (ESR), as well as neutron scattering (NMR), have revealed powerful experimental 
techniques to test magnon excitations in quasi-one-dimensional Heisenberg antiferromagnets. In applied magnetic 
field, ESR experiments exhibit thermally activated resonances, depending on the orientation, that could be interpreted 
as transitions between magnon states whose energies are split by magnetic fields and crystal-field anisotropy. The ESR 
power absorption is proportional to the imaginary part of the dynamical susceptibility I(lu) cx w/mx(w). We calculate 
the ESR transition frequencies in two cases, with applied magnetic field h along the DM vector and perpendicular to 
it. The first case is straightforward, the field couples to the doublet ^^ 2 and ESR intensity is non-zero only when 
Ljfi = h, where u> is the microwave frequency. From the point of view of ESR spectra the resonance has zero width. 
The result holds both in presence and not of DM interaction, so that a measure of ESR spectra with applied field 
along D does not permit to disregard the anisotropic from the isotropic interactions. The analysis of ESR spectra for 
the field perpendicular to DM vector is more involved. We will begin with the case in which Dzialoshinskii-Moriya 
interactions exist only along the rungs. The field dependence of the dispersion relation is no longer determined by 
symmetry considerations alone. The Hamiltonian is given by: 



H 



1 a 



4 
a=l 



E 

>=L,R 



(48) 



Compared to the case without DM interaction, the field couples to both singlets £ 3 and £ 4 so that also singlet modes 
participate to ESR resonances. To diagonalize the Hamiltonian (^) we express it in Fourier space as: 



(49) 



k>0 



where: 
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C(fc) = 



/e R (k)\ 
em 
e R {k) 
e L (k) 

\ti(k)J 



= —= / dx((x) 



(50) 



and: 



H(k) = 



t 


uk 


im 


ih cos a 





z/isina 









—im 


—uk 





z/icosa 





ihsma 




— ih cos a 





uk 


im 3 















— ih cos a 


—im 3 


—uk 












— ihsma 











uk 


inn 




\ 





—ih sin a 








—im^ 


—uk 


1 



(51) 



The spectrum of Ti. is obtained by solving the equation det(e — H) — 0. This equation reduces to the form: 



V 3 + a 2 y 2 + aiy + a 







(52) 



where y — e(k) 2 , and whose coefficients are listed in the Appendix In presence of the KSEA interaction we 
have a slight simplification coming from the condition mi — mi. This permits a reduction of the matrix to blocks 
simplifying the diagonalization. Some asymptotic expressions for the spectrum are obtained in the large and small 
h limit and they are explicitly described in the Appendix 0. As discussed in Ref. |34|, the ESR signal will result 
from transitions between the two lowest branches e i (k 1 a,h){i = 1,2). Transitions will occur at values of field and 
wavevector satisfying \ei(k, a, h) — e 2 (k, a,h)\ = ivh, for a fixed DM interaction. In the anisotropic case the difference 
will depend on momentum k, and a causing a resonance broadening. Comparing with the results of ESR analysis of 
the one-dimensional Heisenberg antiferromagnetEj, in presence of a pure DM interaction the three magnon branches 
mix even for a field along one of the symmetry axes due to the easy-plane effect. If DM interaction would not be 
present a branch would be unaffected by the external field, and only the other two would split. ESR experiments with 
field orientation along symmetry axes would help to confirm the picture of a DM induced mixing. As shown in Fig. |^, 
under the application of the magnetic field perpendicular to the DM vector, a remarkable result is that in contrast to 
the case of an isotropic ladder, the magnetic field does not close the gap, thus gapned but incommensurate correlations 
develop. This result is in consistent agreement with neutron diffraction studiesEj on CaCu203, whose structure is 
similar to that of prototype two- leg spin ladder compound SrCu203. The results show that the magnetic structure 
is incommensurate in the direction of the frustrated interchain interaction in presence of an applied magnetic field 
along the a-axes. An evaluation of the DM interaction reveals that DM coupling in this system may be as large as 
several meV. In presence of a KSEA interaction, the numerical results show a closure of the gap upon application 
of the magnetic field. The KSEA interaction makes the axis parallel to D equivalent to the one perpendicular to it, 
restoring the SU(2) symmetry. This result is shown in Fig. 0. The most general orientation of the magnetic field in 
the plane perpendicular to D will be dealt elsewhere. 

We now briefly discuss the case in which we have Dzialoshinskii-Moriya interaction both along the rung and along 
the chain direction. Then, we have to perform a gauge transformation (Q). As a result, the coupling to the field 
parallel to x direction becomes: 



Afield = -hy^(cos(8n)S% - sm(Sn)Sl) 



(53) 



Using the expressions ( p9| ) of the uniform component of the spin density, we obtain 

'5x\ 



-ff field = — ih J dx(cosa$ — sina^ 



cos 



(54) 



Thus, the magnetic field is now acting as a periodic potential on the Majorana fermions and band structure in the 
spectruni|-should be expected. The discussion of the effects of such band structure on ESR response will be discussed 
elsewhereE3. 







VI. CONCLUSIONS 



We have presented a Majorana fermion description of spin 1/2 ladders in presence of anisotropic superexchangc 
spin interactions along the rungs which allows for the calculation of static and dynamic magnetic properties. We 
have shown that in presence of a pure DM interaction, the systems always has a spectral gap and the lower lying 
excitations are Majorana doublets, instead of triplets for the regular ladder, indicating that the DM interaction alone 
is acting as an easy-plane effect. We have computed both the asymptotic behavior and temperature dependence of the 
anisotropic spin susceptibility helpful for NMR experiments. Our analysis reveals that the easy-plane effect induced 
by a pure DM interaction gives shorter correlation lengths along the DM vector, compared to the orthogonal direction. 
The symmetric counterpart of the DM interaction, i.e. the KSEA interaction, restores the SU(2) symmetry and the 
same correlation length is found in all direction. We have demonstrated that the interplay of the DM interaction 
and KSEA interaction results in very interesting behaviors of the static susceptibility and ESR spectra that are of 
direct experimental relevance. In presence of a pure DM interaction gapped but incommensurate correlations develop 
while the effect of a superimposed KSEA interaction is the closure of the gap when an external magnetic field is 
applied along one of the symmetry axes. The overall analysis reveals that there is a net way of determining the DM 
interaction and the interplay of the DM and KSEA interactions in experiments. It would be interesting to extend the 
present analysis to the strong-coupling limit, or analyse the effect of a canted antiferromagnetism over the DM and 
KSEA interactions. 
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APPENDIX A: UNIFORM COMPONENTS OF THE CORRELATION FUNCTIONS 



According to bosonization, we can decompose the spin density into a uniform and a staggered component as: 

S P = J;?, P + Jl, p + e l ™/ a n p (Al) 

We have discussed the staggered component in the previous sections. In the present section, we derive the expression 
of the uniform .component of the spin density as a function of the Majorana fermion operators. Using the definition 
of the currentsEj and the transformation (^), we obtain the following expressions for the uniform component of the 
spin density: 

{Jm +Jr2) 1 = -i(cosa^| - sina^f^) 
(Jill + JR2) 2 = -i(cosa^R + sina^Cfi) 
(Jill + J R 2? = 

(A2) 

And: 



(Jli + JL2) 1 = -iicosa&l ~ sinaClCl) 
(Jli + Jl2? = -i(cosaClei + sma&£) 
(Jli + JL2? = 

(A3) 



For the differences of the currents, we obtain: 

(Jill - JR2) 1 = -i(cosa^R + sina^^l) 
(Jiii - Jii2) 2 = -i(sbxa^% + cosa^Cfl) 

(Jiii - Jr2? = i&& (A4) 



(Jli - JL2) 1 = -i(cosa£i£ + sma^£) 
(Jli - Jl2? = -i(sina^ 3 L + cosa&ti) 

(Jli- Jl 2 ) 3 ^l£l- (A5) 



From the previous expressions we evaluate the slowly varying component of the Matsubara spin-spin correlation 
functions which are useful for neutron scattering experiments and NMR relaxation rate as reported in the main text. 



APPENDIX B: DETAILS ON THE CALCULATION OF THE NMR RELAXATION RATE AND 

MAGNETIC SUSCEPTIBILITY 



To obtain the NMR relaxation rate or the magnetic susceptibility we have to calculate the following quantity 



T 



dq , Imr a ^(q,w) 
— hm 

2-7T w— >0 LU 



(Bl) 



where T al3 (q, u>) is defined in Eq. (|3^). 
We use the following decompositions: 



GR R {k,iui n ) = -- 



uk 



1 



(k) 



[ 1 



uk 



1 



(k) 

1 1 



e a (k) 



e a (k) 



(B2) 
(B3) 



To obtain after performing the Matsubara sum: 
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T a p(q,u n ) =T™p(q,w n ) +T^(q,uj n ), 



dk 

4tt 



2 k{k — q) + m a rap 
f- a {k)ep(k - q) 



(1 - njr(€a(fc)) - n F {ep(k - q))) 



iuj n + e a (k) +e@(k- q) iu n - e a (k) - ep{k - q) 



iw„ + ep(h — q)- e a (fc) iu n + e a (k) - e fJ {k - q) 



(B4) 



(B5) 



(B6) 



1. NMR rate 



From Equation (B4), it is clear that only T ph (q, ui n ) can make a non-zero contribution to ImT a p for ui — > 0. Making 
the analytic continuation iu> n — > + iO and considering the imaginary part, we obtain: 

,2; 



T 



u 2 k{k~q) +m a m p \ 
1 + awl \ n F e a (fc)) - n F {ep{k - q))) 



x [6(u>- ep(k -q) + e a {k)) - 5 (cj - e a (fe) + ^(fe - g))] 



(B7) 



Performing the integration over q, and making the variable change k — q — ► fc', we obtain the following expression 



for the contribution of T„ n to the NMR rate 
dk \ dk 



T 
8^ 



- a/3 
,2 1. 1./ 



, ( u 2 kk' + m a m \ , S(uj - e (k') + e a (k)) - S(uj - e a (fc) + e p {k')) 
\ l+ e a {k)cf,{k>) J(^M*))-M^(*))) (B8) 

The term containing fcfc' vanishes by symmetry. Since the rest is invariant under k — > — /c or fc' — > — k' , we can reduce 
the integration on (fc, fc') to [0, oo] x [0, oo] and multiply by a factor of 4. It is now convenient to introduce the variable 
change: e = e a (k), e' = ep(k') to find that the contribution to the NMR rate reads: 



T 
2^ 



de 



de' 



^e 2 -m 2 a J mfi ^J(e') 2 - m 



(ee +m a m l 3)(n F (e)-n F {e)) 



For m a > mp and to > 0, this expression is rewrittei 

T 1 f°° 

lim-lmr5( q) .) = -/ de 



30. 



m a m fj 



=cosh 2 ( — ) 

^(e 2 ~m 2 a )(e 2 -m 2 ) 2T 



(B9) 



(BIO) 



For m a — mp, some extra care is needed as the NMR relaxation rate contains a logarithmic divergence for to — ► 0. 
This time, the contribution to the NMR relaxation rate is: 



T 1 f°° 

-Imr ph (q,cj) = — / de 



: cosh 



V2T/ 



Since the divergence comes from e ~ m, to estimate the singular part we can approximate the integral by: 

e~^ T de 



n Jm a m a )(e + w - mj 

Provided m^> T. With the variable change: 



e = m + — (cosht — 1) 



The integral is rewritten: 



771 
7T 



2 e -K-w/2)/T 



die 2T 



— tt^t cosh t 



m 

IT 



£ e -m/T 



0(1) 



(Bll) 



(B12) 



(B13) 



(B14) 



Expansion of the Bessel function for uj —> leads to the formula (3.9) of Ref. 33 
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2. Susceptibility 



The calculation of the susceptibility for m a = mp is straightforward. To obtain the susceptibility for m a 7^ mp, 
we must note that there are two contributions, one from T pp , the other from T ph . Interestingly, the contribution for 
T pp does not vanish for T = 0. Below, we give a detailed calculation of F pp . Starting from Eq. (B5) for ui n = 0, we 
obtain: 



T pp (q^0,io n = 0) 



dk 
2^ 



/ _ (uk) 2 + m a m \ 1 - n F (e a (k)) - n F (e (k) 
\ e a (k)e (k) J e a (k) + e((3(k)) 



Multiplying numerator and denominator by e Q — ep, and taking T — 0, we obtain: 

1 



F*(g->0,£«; n = 0) = 



dk 

2^ 



n ^ | (uk) 2 + m a mp , , (uk) 2 + m a mp 

£u\K) H 7TT £ M J 



To calculate the integral (JB16J) , we consider 

r-K 



e a (fc) 



(ufc) 2 + m a mp 



where K is a cut-off. Using the variable change: uk = m a cosh#, we find: 



K 



dk 
2^ 



e a (k) 



(uk) 2 + m a m/3 



1 

7TU 



uKsJ(uK) 2 + m 2 + m a mp In 



/2uK\ 



With this result, we find that: 



(B15) 



(B16) 



(B17) 



(B18) 



0,w„ = 0) = 



7 — ^ 2T nm 



uK(m 2 a — mp) 



y/(uK) 2 + m 2 a + J(uK) 2 + m 2 



m a mp In 



(B19) 



i.e. 



r( g = o,^ = o) = — 

iru 



2 mi-m^ \m a J 



(B20) 



APPENDIX C: ENERGY LEVELS OF THE DM LADDER IN A MAGNETIC FIELD 

The secular equation that gives the eigenvalues of the Hamiltonian ( |l8|) is of the form: 

y 3 + a 2 y 2 + a x y + a Q = (CI) 
where y ~ t(k) 2 . We put e,(fc) = k 2 + mf. The coefficients ai are given by 



a 2 = -(e\ + e 2 2 + el + 2h 2 ) 
at = h 4 + 2h 2 A(m, a) + ejej + ejej + ejej 

a = -h 4 B(k, m', a) + 2h 2 C(k, m, a) - t\t\t\. (C2) 



A(m, a) = (to 2 cos 2 a — 77717774 sin 2 a + m 2 sin 2 a — m\m?, cos 2 a) 
B(k, m', a) = (to 2 cos 4 a + 277737774 sin 2 a cos 2 a + to 2 sin 4 a + fc 2 ) 
C(A:, m, a) = (mirrism 2 cos 2 a + k 2 m 2 cos 2 a -I- TO1TO3TO4 sin 2 a 

+ k 2 m\m^ sin 2 a + fc 2 m 2 sin 2 a + k 2 m,\mz cos 2 a + fc 4 ), (C3) 

where m = (7711,7713,7714) and m' = (7773,7774). For some specific values of h the solutions of (CI) are plotted in 



Figs.tM. For small h or large h we derive the following asymptotic expressions for the solutions: 
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1. Behavior in strong magnetic field 



For strong magnetic field, it is convenient to search for solutions of the form: y(k) — h 2 u(k). Doing such replacement, 
we obtain the equation for u in the form: 

, ^2 4+4 + 4 2 44 + 44 + 44 \ (B G e 2 e 2 e 2 \ 

"("-i) -^-^ u w ) u -{i^ + ^-n^)- (C4) 

Thus, one can search for solutions of the form: 

u = 1 + E/h + 0{l/h 2 ) (C5) 

u = E/h 2 + 0{\/h 2 ) (C6) 



In particular, we find the solution u = B/h 2 , which leads to e(k) = \f~B~ which is an effectively non-magnetic mode. 
This non-magnetic mode can be obtained in the numerical calculations for high magnetic field (h> 9). We note that 
for a = we recover the singlet mode of the spin ladder. The other solution leads to an energy dispersion of the form: 

e(k) =±7(fe/2,m,a)±/i, (C7) 



where r y(k/2, m) = y k 2 + (rn + 777,3 cos 2 a + sin 2 a) 2 /4. 

We note that for a — and m = m 3 we recover the magnetic modes of the spin ladder. So, in high fields, the 
Dzialoshinskii Moriya ladder appears to behave similarly to a regular spin ladder albeit with renormalizcd gaps. 



2. Behavior in small magnetic field 



At small h, we can expand equation (CI) in powers of h and find a solution valid to first order in h . The equation 
expanded to first order in h 2 reads: 

Y[ (e 2 -e 2 )-2h 2 e 4 + 2h 2 A(m,a)e 2 + 2h 2 C(m,k,a)=0. (C8) 

i=l,2,3 

One can search for solutions of the form: e 2 (k) — ef(k) + fiiti 2 . We find: 

_ [C(fc, m, a) + Ajm, aK(0) 2 - e t (0) 4 ] 

^ ei( o)n^( £l (o) 2 - ej (o) 2 ) ' 1 ' 

From the knowledge of the spectrum we evaluate the total magnetization through the following relation: 
M-X- ST Se *W_ V- V- h [~ei(0) + A 2 (m,a)e 2 (0) + C(k 1 m,a)] 3 

i=le,(fc)<0 i=le,(fc)<0 V ' " l A ' 1 X " 

We can perform analytically the integral and the result gives a finite magnetic susceptibility in zero magnetic field: 
^ [m 2 m 2 ln(^) + (m 2 + m 2 m 2 - to 2 ) ln(^) + (m§ + m|ml - mj) ln(^)] 

Xl'JJ — ( 1 2\( 1 2\ I 2 2*\ ' \ J 

Thus, at zero temperature, there is a finite susceptibility in the direction perpendicular to the DM vector, but zero 
susceptibility parallel to the DM vector. 

In presence of the DM+KSEA interaction the equation ( p8p still holds with the prescription e± = £2 and m = m.3, 
thus we search for solutions of the form: e 2 (k) — e 2 (k) + a.\h + a>2h 2 and e 2 (k) = ef(fc) + a^h 2 . We find: 



[C(k, m, a) + A(m, a)ei(0) 2 - £l (0) 4 ] 

ai -V ( £l (o) 2 - £3 (o) 2 ) ;a2 ~° 



a2 = V ( £l (o) 2 - £3 (o) 2 ) ;ai = 

[C(fc,m,Q) + A(m,a) e3 (0) 2 -e3(0) 4 ] 
£3 (0)(e 3 (0) 2 - e i(0) 2 ) 2 ■ 
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The first solution gives zero contribution to the net magnetization, thus we find: 



m= ^ hem^i (ci( o) 2 _ C 3 ( o)2) — 7m (|oRf + oih) - (C13) 



ei,3(fe)<0 



The integration of such expression gives the zero temperature susceptibility in presence of the DM+KSEA interac- 
tion. 
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D=0 
D=0.8J ± 




T 

FIG. 1. A plot of specific heat versus temperature for a spin ladder with and without DM interaction. A remarkable feature 
is that at low temperature, the specific heat of the ladder with DM interaction is lower while it is higher at high temperature. 




T 



FIG. 2. A plot of spin susceptibility along the z-axes versus temperature for a spin ladder with DM interaction and DM plus 
KSEA interaction, for D/J± = 0.8. The KSEA interaction restores the SU(2) symmetry but a larger spectral gap appears in 
the low-T limit. 



1G 




FIG. 3. A plot of the dispersion of the spin ladder with DM interaction for fixed magnetic field h = 3, u = 1; m is the 
Majorana doublet mass, mz,4, are the masses of the singlets . We notice that there is no closure of any gap upon application of 
a quite strong field. From the dispersion, we should expect gapped but incommensurate correlations to develop in the system. 



m=m 3 =l 

m 4 =3 

a=7t/6 




FIG. 4. A plot of the dispersion of the spin ladder with DM+KSEA interaction for fixed magnetic field h — 1.5, u = 1. We 
notice a closure of the gap upon application of the magnetic field. Thus, the KSEA interaction makes the x axis equivalent to 
the z axis. 
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